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Suppression of resistive hose instability in a relativistic 
electron— positron flow 
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ABSTRACT 

This paper presents the effects of electron-positron pair production on the linear 
growth of the resistive hose instability of a filamentary beam that could lead to snake- 
like distortion. For both the rectangular radial density profile and the diffuse profile 
reflecting the Bennett-type equilibrium for a self-collimating flow, the modified eigen- 
value equations are derived from a Vlasov-Maxwell equation. While for the simple 
rectangular profile, current perturbation is localized at the sharp radial edge, for the 
realistic Bennett profile with an obscure edge, it is non-locally distributed over the 
entire beam, removing catastrophic wave-particle resonance. The pair production ef- 
fects likely decrease the betatron frequency, and expand the beam radius to increase 
the resistive decay time of the perturbed current; these also lead to a reduction of the 
growth rate. It is shown that, for the Bennett profile case, the characteristic growth 
distance for a preferential mode can exceed the observational length-scale of astro- 
physical jets. This might provide the key to the problem of the stabilized transport of 
the astrophysical jets including extragalactic jets up to Mpc (~ 3 x 10 24 cm) scales. 
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1 INTRODUCTION 

Supra-parsec-scale transpo rt of astrophysical jets has been a 
puzzle for a long time (e.g. lBridle fc Perlevlll984r ). and sub- 
stantial efforts have been devoted, both theoretically and 
observationall y, to sol ve their extremely stable feature (for 
a review, see lHarded 12003). Accordin g to the first princi- 
ple analysis (|Appl fc Camenzi nd 1992), the outflows driven 
by a central engine are expected to carry a huge current 
that could reach 10 17 A and more (for the jet of 3C273, 
see, for example. IConwav et al.lfl993h . However, it is known 
that reflecting the interplay between the beam and the self- 
gene rated azimuth al magnetic field, the net current is inhib- 
ited (|Alfvenlll939l ). The upper li mit associated with a unit 
current is given by (lHondall2000l ) 



28/?r kA, 



(1) 



where F — (1 — /3 2 )~ 1//2 is the Lorentz factor of an elec- 
tron (or positron) flow. Apparently, here there is the se- 
rious difficulty that the current of jets greatly exceeds 
7o for conceivable F- values. One possible idea to over- 
come this is to make the cluster of current filaments, 
each having a current belo w 7o, carry such a huge cur- 
rent (|Honda fc Hondall2002l ). This situation might be ac- 
complished via the current filam entation insta bility (CFI) 
including the Weibel mode (cf. iHondal |2004| . and refer- 



ences therein). Indeed, the perfect conductivity character- 
istic of plasma calls the total return current equal to the 
outgoing current, likely arranging a pattern of the coun- 
terstreaming currents, which could serve as a free-energy 
source to amplify the perturbed magnetic fields transverse 
to the currents. In an environment with body-wave per- 
turbations unstable for such an electromagnetic mode, a 
two-stream i nstability for lo ngitudinal electrostatic pertur- 
bations (e.g. iBunemanll 1958T ) could simultaneously develop, 
although this is considered not to play an es sential role in es- 
tablis hing the clustered current system (e.g.lNishikaw a~et al.l 
2003). Also, highly self-organized filaments are certainly 
confirmed i n many radio sources, for example the Galac- 
tic Centre ilYusef-Zadeh et al.l fl984l ; lYusef-Zadeh fc Morris] 
ll987l ; |Yusef-Zadeh et alj|2004 k CygA (iPerley et all Il984l ) , 
M87 lOwen et alJll989h. 3C353 (ISwain et al.lll99ct). 3C 27 3 
(|Lobanov fc Zensusll200ll ) and 3C438 l|Treichel et al.ll200ll ). 
Moreover, it ha s recently been suggest e d tha t the filamen- 
tary jet model (|Honda fc Hondal |2004 l2005f ) could repro- 
duce th e synchrotron X-ray spectra observed in extragalac- 
tic jets jFleishmaj2006l ; iHonda fc Hondall2007l ) . 

Long-scale jets have narrow opening angles less than 
10°, although in close proximity to the central engine the 
an gles tend to signifi cantly spread (e.g. ~ 60° for the M87 
iet: |Junor et all 1999k In addition, there is a consensus that 
the internal pressure of the jets is higher than the exter- 
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nal pre ssure: for example 4C 32.69 llPotash fc Wardlelll980l) . 
CygA i|Perlev et all 1 1984 ) and M87 l|Owen et al.lll989l ) It 
is implied that the jets are self-collimating. In a promising 
case where the outgoing currents are collimated by the self- 
generated magnetic fields, the magnetic pressure exerts to 
preferent ially evacuate the p lasma return currents radially 
outward l|Honda et al.ll2000al lbh. The resulting radial expan- 
sion might be pronounced for the plasma around the enve- 
lope of the filament cluster, because of the lower external 
pressure. As shown by the observational facts, a portion of 
the evacuated flows can constitute a halo surrounding jet 
(e.g. 1803+784; iGabuzda fc Chernetskiil|2003l; iBritzen et all 
120051 ). filamentar y back streams (3C 84; lAsada et alj|200rj) 
a coc oon (Mr k 3. icapetti et al.|[l999l; 3C 273. iBahcall etall 
1 19951 ; CygA, ICarilli fc Barthell Il99fj ). and so on. Impor- 
tantly, as a result of the continuous evacuation, the vacuum 
region s could appear to m ediate the beam and return cur- 
rents l|Honda et al.|[2000al lbl): the velocity shear- free config- 
uration is likely established. Then, the surface mode insta- 
bility driven by the excess kinetic energy of the shear flows 
(i.e. the Kelvin-Helmholtz instability i n a fluid context; e.g . 
iFrank et al]|l996l ; lHanasz fc Solll996l ; iMalagoli et"allll996l ) 
is not crucial. However, if the fractional return current re- 
mains in the beam core, the instability driven by the coun- 
terstreaming currents, such as the hose instability for kink- 
type perturbations, can grow. The eigenfunction is peaked 
near the surface, involving, for example, for a simple beam 
with rectangular radial density profile, the collective reso- 
nance with the beam electrons undergoing betatron oscilla- 
tion (cf. Section 5.1). In this sense, the hose instability might 
also be classified into the surface mode instability. The bulk 
mode instabilities le ading to sausage- or hollowing-like de- 
form ation of beam (|Uhm fc Lampel Il98ll ; iJovce fc Lampei 
1 19831 ) are ignored here, to isolate the axially asymmetric 
modulation. 

The beam kinetic energy tends to be, in part, con- 
verted to transverse therm al energy via collisionless pro- 
cesses l|Honda et al.ll2000bh , whereupon the pressure equilib- 
rium is hydrodynamically determined by the transverse dy- 
namics (|Hondal2000l ; iHonda et al1l2000al ) . However, the con- 
ventional magnetohy drodynamic (MHD) description (e.g. 
IChand rasckhar Il98ll ) is insufficient for the stability anal- 
ysis, as it averages out the betatron oscillation of beam par- 
ticles, merely remaining an averaged flow velocity. If any, the 
cold return plasma that comprises gyrating particles can be 
modelled as an MHD fluid (including the negligible effects 
of particle inertia). However, note that, in realistic situa- 
tions, the energy flux of beams is expected to be larger than 
that of the return flows, such that the beam dynamics dom- 
inantly influences the evolution of the entire system. As a 
consequence, for a collisionless, isolating beam (surrounded 
by the vacuum) the Vlasov equation is found to be ade- 
quate for resolving the wave-beam particle interaction. In 
a fluid-like way, the beam column appears to be stable for 
the velocity shear-free configuration. It is also noted that 
macroscopic dynamics of the cluster of filaments could de- 
viate from that of a uniformly filled cylinder described by 
the fluid equation, as the filaments are electromagnetically 
decoupled, retaining the coherence of the betatron orbital 
motion (reflected in the beam structure). According to the 
kinetic simulations of asymmetric counterstreaming currents 
(reflecting the energy flux budget mentioned above) , indeed, 



the mutual coupling of filaments seems likely to be poor in 
the fully developed non-linear phase of the CFI (except for 
a peculiar phase of rapid c oalescence of filamen ts), owing to 
the return current sheath (jHonda et al.ll2000al ). In another 
regime in which filamentary turbulence is well correlated 
all over the system (larger than the amplitude of the be- 
tatron oscillation and the gyroradius), the kinetic analysis 
could provide a microscopic basis for the macroscopic fluid 
approach particularly inv olving an anomalous re sistivity for 
the filamentary medium (|Honda fc Hon da 2005). However, 
one should notice the crucial point that the collective wave- 
particle resonance possibly gives rise to a morphologic catas- 
trophe in the non-linear phase; this is sheer unpredictable 
in the fluid context. 

In this aspect, the recurrence of the pioneering kinetic 
analysis bv lUhm fc Lampel l|l980l ) may be useful for study- 
ing the fundamental transport property, although their ar- 
gument was limited to be of a stability check of the elec- 
tron beam (without positrons) produced in the laboratory. 
The key consequence is that for the self-collimating elec- 
tron beam with a diffuse radial dens ity profile refle cting the 
so-called Bennett pinch equilibrium (|Bennettlll934l ). the res- 
onant divergence of a growth rate is removed, in contrast to 
the aforementioned simple case with a rectangular profile . 
Related to this, the kinetic simulations l|Honda et al.ll20"00al ) 
reproduced t he self-organization of a Bennett-like profile 
l|Hondall2000l ). Thus, the stable (or quasi-stable) transport 
of self-collimating flows is anticipated to be realized in vari- 
ous situations exhibiting similarity. In particular, for an ap- 
plication of their analysis to astrophysical jets, it is very 
impor tant to take account of t he positron abundance effect s 
(e.g. iRoland fc Hermsen1ll995l) - BL La cs ilXie et al.lll995l) 



M87 ^Reynolds et al.lll996l ) and 3C279 (jWardle et al 



1998) 



- whose details have not substantially been investigated so 
far. 

In this paper, the Vlasov-Maxwell analysis of re- 
sistive hose instability , which was briefly introduced in 
IHonda fc Hondal l|2002l ). is expanded to systematically sur- 
vey the possible unstable modes involved in the relativis- 
tic electron-positron flows with a Bennett-type profile and, 
for comparison, the rectangular profile. In the astrophysical 
context, the focus here is on highlighting some noticeable 
features of pair production effects, rather than a parame - 
ter survey such as was performed bv lUhm fc Lampel (l980). 
It is shown that the charge screening effects by positrons 
directly lower the betatron frequency to suppress the lin- 
ear growth of the convective mode. The lowered betatron 
frequency appears to be reflected in the expansion of the fil- 
ament radius, which prolongs the resistive decay time of the 
perturbed current, to suppress the growth of the absolute 
mode; whereupon for the Bennett profile case, the spatial 
growth distance can be safely longer than (or comparable 
to) the observational length-scale of astrophysical jets. For 
this case, for convenience, the scaling of the growth distance 
in the relevant plasma-parameter region is explicitly shown. 

In order to present these details in a straightforward 
way, this paper has been organized as follows. In Section 2, 
an ad hoc model of negatively charged electron-positron 
gas is introduced. For the relativistic flows, in Section 3 the 
equilibria of the momentum distribution functions and self- 
generated magnetic fields are constructed for the cases of the 
rectangular density profile (Section 3.1) and Bennett profile 
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(Section 3.2). Then, in Section 4, the Vlasov-Maxwell equa- 
tions are linearized around the equilibria. Using the pro- 
cedure outlined in Appendix, the eigenvalue equations are 
derived from the linearized equation, to extract the relevant 
eigenmodes in Section 5, for the cases of the rectangular 
(Section 5.1) and Bennett (Section 5.2) profiles. Section 6 is 
devoted to a discussion concerning mode preference in actual 
circumstances, and a summary. 



the range of < (3 ^ 1. Similarly, the effective temperature 
is defined as 

T = f h T b + (1 - / b )T r , (4) 

where Tb, T r = const. For the special case, ft = 1, reflecting 
that the plasma return current is perfectly evacuated outside 
the beam core (Section 1), equations {[3j) and Q reduce to 
P — Ph and T = Tb, respectively. The definitions given by 
equations d2j— (J4j) are substantially used throughout. 



2 A SIMPLE MODEL OF 

ELECTRON-POSITRON GAS 

The extreme astrophysical environments, in which electron- 
positron pairs are created, can be typically found in g alac- 
tic nuclei fe.g. IZdziarski fc Lightmanlll985l; lHevlll200ll ; for 
a laboratory experiment, see IWilks et al.ll2005l ) As a rule, 
the charge conservation law requires that the total electron 
charge is compensated by the pair-created positrons and dis- 
charged ions, such that n c = n s + {Z*) rii, where n c , n s 
and Hi are the total electron, positron and ion densities, re- 
spectively, and (Z*) is the average of the charge state of a 
composite. Hereafter, the superscript '"' indicates the quan- 
tity for 'anti'-electron (i.e. positron). For simplicity, the ra- 
tion of ric/n c = /pair is supposed to be constant, where 
^ /pair < 1. Then the definition 



A/pair — 1 /pa 



(2) 



is introduced, and the reciprocal, (A/ pa j r ) , is referred to 
as the pair production rate; these take the values in the 
ranges of < A/ pa ir ^ 1 and (A/pair) -1 ^ 1, respectively. 
Because n c > ric, the electron-positron gas may be regarded 
as negatively charged. 

As discussed in Section 1, in the short-circuit reflect- 
ing the characteristic of perfect conductivity in plasma, the 
return currents are allowed to flow outside the beams, an d 
even the jet (for modeling, cf. lBenfordfl97Sl ; lAlfvedfl98ll ). 
For the present purpose, the fractional beam density, de- 
fined as the ratio of the beam electron (positron) density 
to n c (rio), is assumed to be constant [i.e. /b = rib/n e (= 
rib/rie) = const], whereby the fractional return flow den- 
sity is given by 1 — /b = n r /ro c (= nj/ric) = const. Then 
we have rig = /pairrib and rif = /pair^r- Hereafter, the 
subscripts 'b' and 'r' indicate the quantities for the beam 
and return flow, respectively. It is also mentioned that, for 
(A/pair) -1 > 1840 (Z*) I (A*) ~ 10 3 and < 10 3 , where 
{A*} is the average of mass number, the electron-positron 
fluid and ion rest frames, respectively, might be properly re- 
ferred to as the 'jet frame'. In order to exclude the ambiguity 
of the reference, throughout this paper, all physical quanti- 
ties, except for the quantity £1 introduced in equations (1291) 
and (|30l) , are specified in the ion rest frame l|Honda fc Hondal 
I2002T I. 

We consider the situation that the electron-positron gas 
flowing in the z-direction has the speed of Vb and the gas 
flowing in the opposite direction has the speed of \v T \ (<§C Vb), 
and these have the thermal energy of Tb and T r , respectively. 
Then, it is convenient to define the effective parameter of 



p = f h p b - (1 - f h )\pr\. 



(3) 



Here, Pb ~ rtb/c = const and |/3 r | = \v v \/c = const, where c 
is the speed of light. Note that equation @ takes a value in 



3 EQUILIBRIA OF THE RADIALLY 
CONFINED RELATIVISTIC FLOWS 

In what follows, the zeroth-order distribution functions of 
beam electrons are assigned, as they provide, within the 
present framework, sufficient information to characterize the 
equilibrium prop erties of the negatively charged fluid (cf. 
IWilks et alJfeOOa ). First, we consider a simple case of a col- 
limated flow with thermal spread and sharp radial edge (Sec- 
tion 3.1), and secondly, another case of a self-collimating flow 
with thermal spread and diffuse density profile (Section 3.2). 



3.1 The case of rectangular radial density profile 

Using the similar procedure by iDavidsonl l|l990l ). we begin 
with the equilibrium distribution function for the beam elec- 
trons with the momentum of p = (p r ,pg,p z ): 



/b°(p,r) 



rib 



27T7 b m, 



■F h (H)S (p z - 7 b m e /3 b c) 



(5) 



Here, Fb(H) = 5 \ H — (jb — 1) m c c 2 ] , 8 indicates the Dirac 
delta function, H = {jp\ + pj) / (2j h m c ) + (j h — l)m c c 2 + 
e(3bA s z (r) is the Hamiltonian, e is the elementary charge, 
m c is the electron rest mass, 7b = (1 — P^)' 1 ^ 2 , and 7b 
(> 7b) includes the thermal component. In the cylindrically 
symmetric case, the z-component of the vector potential, A%, 
conforms to the ^-component of the magnetic field: Bg(r) = 
— dA z (r)/dr, where the superscript 's' indicates the 'self- 
generated quantities. Equation ([5} is rewritten as 



fb(l>,r) = 



rib 



27T7bm e 



27b m< 



+ 4>b(r) - (7b - 7b) m c c 

x8 (p z - 7 b m c /3 b c) , (6) 

where p\ = p\ + pg, and ipb(r) = ePbA B z (r) represents the 
effective potential that acts on beam electrons. 

For the definition of the beam electron density, rib(r) = 
J d 3 pfb(p,r) = ty / °°d(pi) f™ ca dp z f h ) (pj_,p z ,r), we have 



n b (r) = h h / dUS [U + V>b(r) - {■% - 7b) i 



(7) 



where U = p 2 L /(2jbm e )- The integral is equal to unity 
for ipb(r) < (7b — Jb)m c c 2 and equal to zero for ifjb{r) > 
(7b — 7b) m c c 2 . Therefore, the density profile has the simple 
rectangular form 



rib 



fib = const for ^ r < rb 







for r > Tb 



(8) 



whereby ng = ng — const for ^ r < rb and ng = for r > 
rb- The beam filament radius, rb (smaller than the radius of 
the jet; Section 1), can be self-consistently determined by 
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tp{r h ) = (7b - 7b) m e c 2 . 



0) 



Note here that the spatial profile of the cylindrical magnetic 
field can be expressed as 



Bl{r) = 



-2-7re/3A/pai r n e )" 



for ^ r < rb 



-27re/3A/ pa i r n e (r b /r) for r > r b 



(10) 



where n c = hi,/ fb, an d the corresponding effective potential 
is ipb(r) — 7re 2 /3b/3A/ pa i r n e r 2 (for < r < r b ). Therefore, 
equation © yields 



2 2c 7b - 7b 

r b = — 



(11) 



where u}p(= const) defines the betatron frequency, whose 
square is given by 



.2 27re /3 b /3A/p air n c 
w fl — 



7b«1c 



(12) 



It is found that the larger pair production rate (i.e. the 
smaller value of A/ pa i r ) leads to smaller uip and larger rb. 



3.2 The Bennett-type equilibrium case with a 
diffuse radial density profile 

We consider the radial force balance between the gradient 
of static pressure stemming from the thermal spread of the 
beam, (d/dr) (1 + / pa i r ) n e (r)T, and the J x B contraction 
force. This might reflect a more realistic equilibrium of rel- 
ativistic elec tron-positron flow. When the electrostatic field 
is negligible dHonda et al.ll2000al). th e force balance equation 



written as jHonda fc Hondai 20021 



T(l + / pair ) 



dn c (r 
dr 



P A/ pair n c (r) 



x / dr'r'n c (r'), 
Jo 



(13) 



where we use the expression of the cylindrical magnetic field 
self-generated as a result of the net current of an electron- 
positron flow: 



B B e (r) 



4ne 



/3A/ pa i r / Ar'r'n c (r'). 



(14) 



The spatial profile of the electron number density can be 
self-consistently determined, to have the form of 



n e (r) 



n h (r) 



lib 



1 



(l + r7rg) 



2^2 ' 



(15) 



where h c = n c (r = 0) = rib/ fb, and the characteristic radial 
size of a beam filament is given by 



2 

'b 



2(l + / pair )T 
^ 2 /? 2 A/ 2 ri c 



(16) 



In contrast to equation ([S]) , the density profile given by equa- 
tion (115[) has no sharp radial edge. This corresponds to the 
modified Bennett pinch equilibrium of a self-collimating rel- 
ativistic electron-positron flow, in which spatial profile of 
the magnetic field can be expressed as 



-27re/3A/ pair n e — — . 

1 + r l r b 



(17) 



Note that B B g (r) takes the peak value at r = rb, where the 
magnetic energy density, \B s e (r = rb)| 2 /(87r), is in the level 
of ~ (1 + / pa ir) h B T, reflecting pressure balance. 

Let us consider the equilibrium beam distribution func- 
tion of the form of 



/b°(p,0 



»b 



27T7bm e 



Pi 
27b m, 



+ (7b — 1) m e c 



+e{3 b A B z (r)] 8 (p z - 7 b m e /3 b c) 



(18) 



Recalling the definitions of U and ipb(r), the beam electron 
density is then given by 

/>oo 

n b (r) =h b dUF b [U + ( 7b - 1) m c c 2 + ^ b (r)] . (19) 
Jo 

Imposing the boundary condition of A%(r = 0) = 
0, we integrate equation {XT}, to obtain nb(r)/rib = 
exp [-2A%(r)/ (7re/3A/ pair n c r^)] , which is recast to 



exp 



n b (r) 
n b 

Here, Vb(r) 
dij) b {r 



/3A/ pair 



;^b(r) 



Pb (1 + /pair) T 

ePbA B z (r)] conforms to the equation of 



dr 



7 b m e ra) /3 (r 



(20) 



(21) 



where uip defines the betatron frequency, whose square is 



1 + r 2 /r 



(22) 



Invoking the density inversion theorem for equation (|19[) . 
namely, F b (H) = -n~ 1 [dn b /dil> b ]^= H _ ( _ 1)m , e 2 wher e 
H = U + (7b - l)m c c 2 + Vb (Section 3.1; lDavidson|[l990l ). 
we explicitly obtain the expression of F b (H), that is, 



Fb(H) 



/?A/ pa 



: exp 



/3A/ pa 



Pb (1 + /pair) T ' \ P b (1 + /pair) T 

x[H- ( 7 b - l)m c c 2 ] } . (23) 

This constitutes a pseudo distribution function of thermal 
equilibrium for equation (|18[) . 

If we take account of no direct contribution of re- 
turn current component to the thermal pressure gradi- 
ent and J x B force (in the radial force balance), in 
equations (I13|l and ()16p the replacements of T — > /bTb 
and p 2 — > fbPbP are valid, whereby in equation (|23[) , 
/3A/ pair /[/3b (1 + /pair) T] -> A/ pair /[(l + / pa i r ) T b ]. In this 
regime, for the special case of no pair production (i.e. 
/pair = and A/ pa ir = 1), equation (f23|) reduces to 



T, 



exp 



H - (7b - l)m c c 2 



(24) 



which coincides with equation (9.239) in iDavidsonl |l990). 



4 PERTURBED VLASOV-MAXWELL 
EQUATIONS 

Now we superimpose the fluctuations on the aforementioned 
equilibria of f b and A% , such that the total distribution func- 
tion and vector potential can be linearized as f b = f b + 5f b 
and A = A\i + SA, respectively. Assuming a resistive re- 
sponse of return currents, the fluctuating total current can 
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be phenomenologically expressed as <5J = A/ pa i r <5Jb + crSE, 
where <5Jb denotes the fluctuating beam electron current, 
a is the electrical conductivity and <5E = —(l/c)(d5A/dt) 
is the fluctuating electric field (without electrostatic com- 
ponent; cf. Section 3). We also assume a slow spatiotem- 
poral chang e, such that \ui \ <C 47rq, \uj\r^/c <C 1 and 
\k z \r h < 1 (|Uhm fc Lampei Il980l ; Ibavidsonl Il990l ). Then, 
the linearized Ampere-Maxwell equation can be expressed 
as V x <5B = (4n/c)8J, whose z-component is 



ld_ (d_\ , 1 d 2 

r dr \ dr J r 2 d9 2 



47rcr(r) d 



dt 



SA z (r,t) 



4ne 



A/pair / d p«b5/ b (p,r,t) ; 



(25) 



for the cylindrically symmetric case. Here, the terms pro- 
portional to c~ 2 (d 2 /dt 2 )SA z and (d 2 /dz 2 )8A z have been 
neglected as consistent with the above assumptions. In ad- 
dition, the linearized Vlasov equation is described as 



+v h -^- +e/3 h Bg(r) 



dt 



dz 
= e/3 b 



d 

dp r 



Sf h (p,r,i) 



+ 



l dSA z (r,t) 0/g(p,r) 
r 80 dps 

dSA z (r,t) 9/ b °(p,r) 



dr 



8p r 



(26) 



In order to derive the dispersion relation from above set 
of equations (|25[) and (|26|) . we assume the standard Fourier 
expansions of the perturbed components: 



<5/ b (p, r, t) = 5/ b (p, r) exp [i (£9 + k z z - cut)] . 
SA z (r, t) = 5A z (r) exp [i (£$ + k z z~ at)] , 



(27) 
(28) 



where i = \J — 1 is the imaginary unit. In equations and 
(|28|) . it is convenient to transform the set of the independent 
variables (t, z) to (r = t — z/vb, z), so that 

Sf b (p, r, 9, z, r) = 5/ b (p, r) exp [i (£9 - Q,z/v h - ur)\ , (29) 

5A Z (r, 6,z,t) = 8A z (r) exp [i (£9 - Qz/v h - ut)] , (30) 

where Q, = ui — k z Vb is the frequency o f the perturbations 
seen by a beam electron (|Davidsonf l990). Below, we concen- 
trate on an interesting case imposing the kink-type pertur- 
bation with £ = 1 and seek the corresponding eigenmodes. 



5 MODIFIED EIGENVALUE EQUATIONS FOR 
KINK PERTURBATION ON AN 
ELECTRON-POSITRON BEAM 



Using the analysis by lUhm fc Lampei (|l980l ). we derive 
eigenvalue equations from the self-consistent equations (|25p 
and (|26|) : (i) for the case of the rectangular density pro- 
file (Section 5.1) and (ii) for the Bennett profile case (Sec- 
tion 5.2). 

5.1 The rectangular radial density profile case 

According to the derivations outlined in Appendix, we de- 
rive an eigenvalue equation for the kink-type perturbations 
superimposed on the filamentary flow. For the rectangular 
density profile given in equation (|8}, we have drih/dr = 



~~hbS(r — rb). Then, the master equation (|A9f) provides the 
following eigenvalue equation: 

rid_ d_ 

\r dr dr 
8A Z 



1 47rkj(j\ . , 
-i, + — — I 6 A 



J SA z (r) 



\ ~ 2 
)^pb 



/?b/7b 



(n 2 



S(r- 



(31) 



47re 2 ribA/p a i r /m c = const 



where the definition of i2i 2 b 
has been introduced. It is evident, from the right-hand side 
(RHS) of equation (|31[) . that the perturbed axial current is 
localized at the beam surface (r = rb), analogous to the 
Kelvin-Helmholtz-type instability for neutral flows with lo- 
cal velocity shear (e.g. IChandrasekharl Il98"lh . Apparently, 
eq uation (|31l) has t he same form as equation (9.256) given 
m iDavidsonl (|l990l ). but now the pair production rate is 
included in ui p b and Up. It is expected that in astronomi- 
cal environments with large pair production rates, both a) p b 
and tip are significantly reduced, to lower the intensity of 
the perturbed surface current and resonant frequency of Q, 
respectively. When we impose the boundary condition of 
SA z (r = 0) = = SA z (r = oo) and continuity of SA Z at 
r = rb, the eigenvalue contained in equation f|31f) is found 
to conform to the following dispersion relation: 



UIT d 



(1 - / m ) (1 - X) 



- 1 



(32) 



Here, we introduce the definition of the characteristic resis- 
tive decay time 



Td 



ivar 2 
~2c2~' 



(33) 



and 

,-,2 



Q 2 /w 2 



13' 



and 



also 
-l 



relation of 
= (1 - 



use the 

^b^/(2^ 7b ) = "(1 - / m )-\ where /, 
fb)\(3 r \/(fbf3b), which takes a value in the range of ^ / m < 
1 (cf. equation |3J) . It is noted that equation (|32p is correct 
to leading order in ur d . 



5.1.1 Mode property for real lo and complex Q 

First, we consider the real ui and complex Q case; the spa- 
tially growing/decaying mode may, for convenience, be re- 
ferred to as the convective mode. To avoid confusion, it is 
noted that the growth of perturbations can be viewed in 
the beam frame (not in the wave frame as conventionally 
chosen). In an interesting regime of ujTd — > (i.e. for low 
frequency and/or very high resistivity), the left-hand side 
(LHS) of equation H32|) vanishes, to yield the simple equa- 
tion of (1 — /m)(l — x) — 1, which contains the solution of 
x = — / m /(l — /m). Therefore, for < / m < 1, the value of 
x = k 2 v 2 /uji is real and negative. It follows that k z takes 
purely imaginary values so that k z — ±ifci, reflecting the 
purely growing — ifci and purely decaying +ifci modes (cf. 
equations [57] and I28p . When the current neutral condition 
is perfectly satisfied [i.e. f3 = (/ m = 1)], the value of 
ki = (<2j/3/vb)[/m/(l — fm)] 1 diverges (and thus violates 
the theory), whereas no return current (/ m = 0) trivially 
leads to ki = 0. Clearly, it turns out that the instability 
mechanism in this regime is a result of the repulsion of the 
beam current by the return current w ith / m 7^ (for the 
turd / case, see lUhm fc Lamp"elll980l ). 

We define the growth distance by L = 2n/ki, which 
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scales as 
L = 0.1 



10 



1/2 



1/2 



10" 3 \ 1/2 /l0- 6 cm-^ 1/2 



A /pair 



(34) 



Here equation (fl2]| . the relation of /m/(l — /m) = (/b/3b — 
0)1 ft = (1 — /b)|/9r|//3 and fi r = (1 — /b)n e are used. Impor- 
tantly, as the pair production rate increases, the value of L 
increases, because of the apparent lowering of u)p (cf. equa- 
tion [12]). Even though the pair production effects increase 
L, the predicted L is much shorter than the length-scale 
of major astrophysical jets, particularly, extragalactic jets. 
Thus, the convectively growing mode (for lower uj and/or 
shorter Td) seems not actually to reflect the characteristic of 
the extremely stable transport of the jets. 

5.1.2 Mode property for complex uj and real fi 

Next we examine the complex u and real f2 case; the tem- 
porally growing/decaying mode may be referred to as the 
absolute mode, using conventional terminology. Because of 
real Q, x (= Q, 2 /&%) is real and positive, so that the square 
bracket on the RHS of equation (|32[) is also real. Thereby, 
uj is found to be purely imaginary: oj = iWi. The complex 
wavelength is expressed in the form of k z — ifei + k r , and the 
relation of SI = uj~k z Vb is recalled. Then, we find ki = uji/vb 
and fc r = — fl/vh, which can, respectively, be expressed as 



h = 



l 

VbTd 



and 



(l-/ m )(l-x) 



= -sgn(n)^y^, 

Vb 



(35) 



(36) 



Here, sgn(S7 ^ 0) = ±1. Note that the positive value of ki 
can now be compared to the growth factor for the absolute 
mode, and when x approaches unity below (i.e. x — * 1—), 
h — > +oo. We can confirm that the destabilization involv- 
ing ki > takes place for < x < 1. The fact of the stable 
transport of astrophysical jets over parsec scales seems to be 
inconsistent with the above prediction of the resonant desta- 
bilization. As shown later in Section 5.2.2, the resonance can 
be removed for the Bennett-type profile case. 

5.2 The modified Bennett profile case 



We again recall equation (|A9fl , multiplying through by r5A z 
and integrating over r from to oo. An appropriate choice of 
the trial function of &A Z oc dA% jdr co rresponding to a rigid 
displacement of the magnetic field (cf. lUhm fe Lampe|[l980l . 
for a complete discussion) yields the extended eigenvalue 
equation of the form of 

dA s z (r)] 2 f°° , dA%{r) dJ(r) 



dr ra(r) 



dr 



dr r- 



dr 

2 i 



dr 



P [n*-«a(r)] 



(37) 



Within the present framework (Section 2), the total current, 
J (= Jb + Jf, + J r + Jf), can be expressed as 



J(r) = 



/?A/ pai 

Mb 



-Mr), 



(38) 



where Jb{r) = e/bn e (r)w b . Using equations (fl~5)) , (fT7)l 
and (|38[) . and the variable exchange of £ = (1 + r 2 /r 2 )" 1 
equation (|37l) is cast to 



2 /■! 

mr h (jJ 



2c 2 



d« (1 - 



n 2 + [e^(i-/ b )|AI//3] 



« 2 



•^1 



(39) 



Note that r 2 , and L)p involve the pair production rate, al- 
though the form of equation (1391) itself is similar to equa- 
tion (73) given in llJhm fc Lampd l|l980h . The RHS of equa- 
tion (|39[) can be exactly integrated, to give 

"rd = i (1 - /m)" 1 [G(x; g) + f m ] . (40) 

Here, we introduce the definitions of 



Td = — 

c 

and 
G(x;g) 



"5 



6x 



x + x (1 



hi 



+ iff"" 



(41) 



(42) 



where g = for x < and x > 1; g = sgn(£7) for < x < 
1. Note that when setting the expression of the electrical 
conductivity to cr(£) = <r£ 2 , where a = const, equation (|41|) 
coincides with equation (133[l , 



: Q 



5.2.1 Mode property for real uj and comple 

In parallel to the argument given in Section 5.1.1, we ex- 
amine the real oj and complex fl case, for convenience, re- 
ferred to as the convective mode. In an interesting regime 
of ujT d — > 0, the square bracket of the RHS of equation (|4T7)l 
vanishes, to give the dispersion relation of G(x; g) + f m — 0. 
In particular, the equation of G(x; g — 0) + / m = is found 
to contain the solution of x — x(/m) that is real and, for the 
possible range of / m , negative (as consistent with the choice 
of g = 0). Hence, recalling the expression of x = k z v 2 /oJp, 
the dispersion relation turns out to contain the complex k z , 
which is purely imaginary so that k z = ±ifci, reflecting the 
purely growing — ifci and purely decaying +ifci mode. Note 
here that h = {ujp/vb^-xif™)] 1 ' 2 ■ 

Fig. [T]plots the dimensionless growth factor (— x) 1 ^ 2 as 
a function of / m for 0.01 ^ f m < 1, together with, for com- 
parison, the function of ki{vb/0jp) = [/m/(l — /m)] 1//2 for the 
case of the rectangular radial density profile (Section 5.1.1). 
When the current neutral condition is perfectly satisfied [i.e. 
(3 = (f m = 1)], the value of h = (^/«b)[-x(/m)] 1/2 di- 
verges, whereas no return current (/ m = 0) trivially leads 
to ki — 0. It is, again, found that the repulsive interaction 
between the beam current and the return current is essen- 
tial for driving the instability. As seen in Fig. (TJ the present 
value of ki is smaller than that obtained for the rectangular 
density profile. Namely, the diffuse boundary effects extend 
the growth distance of L — 2n/ki. However, the effects are 
not so significant that the expected value of L is found to 
be much shorter than the length-scale of astrophysical jets 
(cf. equation I34J1 . More importantly, the pair production ef- 
fects significantly lower the value of ujp, thereby decreasing 
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Figure 1. The growth factor of a convective mode, ki(wr^ — > 
0), multiplied by v^/uup as a function of / m , for the rectangular 
radial density profile (equation \E\ dotted curve) and the diffuse 
radial density profile for the modified Bennett pinch equilibrium 
(equation 1151 solid curve). Here, is the resistive decay time 
(equations [33] and [4TJ, v-^ is the beam velocity, ujp is the betatron 
frequency at the beam centre r = (equation I12|l , and / m is the 
ratio of the return/beam current density. 



fci; however, the value of L seems to be still much shorter 
than the observed length of the jets. 



5.2.2 Mode property for complex uj and real SI 

In parallel to Section 5.1.2, we examine the complex uj and 
real SI case, referred to as the absolute mode. Because of 
real S7, x (= Q, 2 /uj^) is real and positive, so that G(x; g = 0) 
is also real. For g 7^ (for < x < 1), both k z and uj are 
complex, which may be expressed as k z = ifci + fc r and uj = 
iuji + uj t . Recalling equations (|40[) and (|42|l , and the relation 
SI = uj — k z Vb, we have ki = uji/vt and fc r — (oj t — Q,)/vb, 
which can, respectively, be expressed as 



fci 



1 G {x)+f n 



VbTd 1 — /m 

where Go{x) = G(x; g = 0), and 



1 

Vb 



6x 2 (1 



(1 - /m) T d 



+ sgn(Q)ujf3\/x 



(43) 



(44) 



On the RHS of equation (|44|l , the first term inside the square 
bracket is ordinarily much smaller than the second term, 
whereupon equation (I44p can be approximated by equa- 
tion (|36|) . The positive value of fci can be compared to the 
growth factor for the absolute mode, but, importantly, equa- 
tion (|43[) does not contain the resonance point of x that 
appeared in the rectangular density profile (Section 5.1.2). 

Fig. [2ja) plots, for the special case of / m = 0, the di- 
mensionless growth factor fci(t>bTd) of equation (|43[) as a 
function of x in the range x ^ 0, together with, for com- 
parison, equation (|35[) for the rectangular profile case. It is 
noted that for f m = 0, the function fci(ubi~d) of equation (|43[) 
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Figure 2. The growth factor of the absolute mode, fci, multi- 
plied by v h r d as a function of Q 2 /uj^, indicating (a) for / m = 0, 
a comparison in between the rectangular profile (equation I35I 
dotted curve) and modified Bennett profile (equation I43I solid 
curve) cases, and (b) the / m -dependence for the latter case. Here, 
SI (= uj — - fc z i>b) is the frequency of perturbations viewed in the 
beam frame. In (a), for the rectangular case the resonant point 
exists at fl 2 /uj% = 1 (vertical thin dashed line). Note that the 
solid curve in (a) and the curve for / m = in (b) correspond 
to the function G(f2 2 / uj ? ; g = 0) (equation 1421 lUhm &c Lampei 
\l98(j) . 



reduces to Gp(x ), recallin g fig. 2 in lUhm fc Lampei (Il980l ) 
and fig. 9.16 in iDavidsonl l|l990l ). As G (x -> 0) = and 
Go(x — > +00) = —1, in equation (|43p we find fa(x —►()) = 
(I'bTd) -1 !/™/!!- fm)] and h(x — > +00) = -(ubTd) -1 ; these 
asymptotic values coincide with those of equation (J35J . As 
seen in Fig. [2^a), in the range < x < 1, the growth fac- 
tor for the rectangular case is always positive, whereas the 
diffuse effects arrange an x-domain in which the stable con- 
dition, fci 0, is satisfied. Note that Go(x — |) =0; in the 
range < x < i, Go (x) > takes the maximum value of 
Go, max = 0.690 at x — 0.272, whereas in the range x > ^, 
Go(x) < takes the minimum value of Go, mm = —3.36 at 
x — 0.949. In particular, for x > 0.635, Go(a;) < —1, so 
that the stable condition, fci < 0, is always satisfied for the 
possible range of / m - 

Fig- Ob) plots, for / m = 0, 0.5 and 0.9, fci(wbTd) (equa- 
tion I43[) as a function of 1. It is clearly found that when 
the return current increases (i.e. with increasing / m ), the 
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maximum (minimum) value of fci(ubTd), which is achieved 
at the common appoint of x\G =G , max ( x \G =G 0tinia ), enor- 
mously increases (decreases). In contrast to the infinite 
divergence at x = 1 involved in equation (|35|l . Go, max 
is found to provide the finite maximum value of fci [i.e. 
fci,max ^ (v b r d y\0.7 + / m )(l - /m)- 1 at x ~ 0.27], where 
equation (|44|) approximately reads fc r ~ — sgn(r2)0.52[2)/3/iib. 
Namely, the growth factor is bounded; physically, this is 
because there is no single frequenc y Q for which the en - 
tire beam is resonant with the wave l|Uhm fc Lampelll980[ ). 
The values of fci !max and k r determine the shortest growth 
distance of wave envelope defined by L m i n = 27r/fci imax and 
the corresponding oscillation wavelength of the carrier wave, 
A osc = 27r/|fc r |, respectively. These are explicitly written as 



Lmin 



Aosc 



7r p b r h a 1 - /m 

C 0.7 + f n 

12wb 

up 



(45) 
(46) 



In the derivation of equation (|45[) . equation (|33|) has been 
employed as the expression of The obtained result sug- 
gests that astrophysical jets can propagate over, at least, the 
distance of ~ L m i n , having the spatial oscillation with the 
wavelength of Aosc (<C L mill ] • 

Let us evaluate equation (|45|l . Concerning the elec- 
trical conductivity, for simplicity we here employ the 
standard Spitzer formula. In the regime of T > m e c 2 
likely for pair plasmas, we have the relation a /as ~ 
[(27r) 1 / 2 /8](m c c 2 /T) 1/2 , where a s ~ T 3 / 2 /(e 2 m c 1/2 In A) 
is the well-known non- r elativistic Sp itzer conductivity 
jBraams fc Karnevl 1 19871 ; iHondal l2003r); that i s. for the 
Coulomb logarithm of In A ~ 10 (e.g. lHuballl994h . the scal- 
ing of a ~ 10 12 (T/10 10 K) s _1 . We also call the relation 
(1 - / m )/(0.7 + / m ) = /3/(1.7/ b /3b - /8). Then, using equa- 
tions (I16p and the above expression of a for T ^> m e c 2 , 
equation (|45l) is found to scale as 



■*>(&) 



1 cm 

n b 



r 



10 10 K 



Mpc, 



(47) 



while for T <C m e c 2 , setting to <r = <xs, yields 



Lmin = 1.1 x 10 



■ ( 10 \ 



1 cm 
n b 



r 



10 s K 



5/2 



pc. (48) 



Here, the enhancement factor is introduced, defined as 



n 



1 - (A/pair/2) 



/3A/ 2 air [l-(/?/1.7/ b /3 b )]' 



(49) 



which takes a value in the range 77 ^ 2/(1.7/b/3b). The strong 
T-dependence of equations (|47[1 and (|48|) arises from that of 
rl oc T, and a oc T (for T > m c c 2 ) or T 3/2 (for T < m e c 2 ). 

In the regimes of A/ pair /2 < 1 and (5/(1.7 f h (5 h ) < 1, 
equation (149 1) can be approximated by ~ l/(/3A/ 2 ail .) 3> 1. 
In particular, for A/ pa i r /2 1, the factor 77 is reciprocal 
to A/ 2 air (i.e. proportional to the square of the pair pro- 
duction rate), reflecting that the radial expansion stemming 
from the charge screening effects by positrons (rb oc A/^ r ) 
significantly increases the value of Td, to extend L m i n . For 
f3 < 0.85/b/3b, r\ increases with decreasing /3, and particu- 
larly for (3 <C 1.7/b/3b, f] oc /3 _1 ; this implies that the radial 
expansion effect (rb oc /3 _1 ) tends to surpass the destabi- 
lization by the return current. Even for the special case of 
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Figure 3. The contour plots of L m in/'7 for the values of 10 — 8 — 
10 6 pc (labelled) in the rib — T parameter space, for the Bennett 
profile case in Fig. [2] (exhibiting the maxima fc; ma x)- Here, L m i n 
is the shortest growth distance (defined by L m ; n = 27r/fc; max ; 
equation I45I I. r\ (> 1) is the enhancement factor (equation 1491 . 
rib is the beam electron number density at r = 0, and T is the 
effective temperature (equation [4} . The shaded area indicates the 
(rib, T)-parameter domain where supra-Mpc-scale transport is 
allowed for 77 ~ 1. 



/b = 1 and P = Ph (fm = 0), equations (|47]l and (|48|) are 
valid, having V = (2.4//3 b )[l - (A/ pair /2)]/A/ 2 air (> 1.2). 

Fig.[3]plots rib as a function of T, for the given parame- 
ter Lmin/77 in the range 10 -8 — 10 6 pc. Here, the generic for- 
mulae of <j is used, which covers the m ild relativistic region 
l|Braams fc Karnevl [l987l ; lHondall2003r ) and In A dependent 
weakly on density and temp erature (e.g . for rib ^ 1 cm -3 
and T < 10 10 K, In A < 30; lHubdfl994r ). The shaded area 
indicates the (rib, T)-parameter domain in which Mpc-scale 
transport of jets is safely achieved even for the marginal 
(most pessimistic) case of 77 ~ 1. Such a parameter domain 
more significantly expands for the larger values of rj owing 
to, for example, the larger pair production rate. As a con- 
sequence, it can be claimed that there is a wide parameter 
domain in which relativistic outflow can safely propagate up 
to 1 Mpc corresponding to the observed longest scale of ex- 
tragalactic jets. Also, for fixed 77, the smaller Lmin allows the 
flow to possess lower T and higher rib, that is, expanding the 
allowable (rib, T)-parameter domain. This property seems to 
be reflected in small-scale jets such as Galactic jets, compris- 
ing lower temperature, electron-ion plasma with T m c c 2 

and / pair = (A/p a i r = 1). 



6 DISCUSSION AND CONCLUSIONS 

As investigated above, absolute growth can be allowed, par- 
ticularly for the case of the diffuse radial density profile re- 
flecting the Bennett pinch-type equilibrium. However, in a 
low- frequency and/or high-resistivity regime, the character- 
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istic distance of the convective growth is, for both the rect- 
angular density and Bennett profiles, much shorter than the 
observed length of astrophysical jets, even if the pair pro- 
duction rate is large enough to lower the betatron frequency 
Cop. Hence, the actual development of the convective mode 
appears to be ruled out. A simple explanation for the rea- 
son why the mode could not crucially evolve is that the 
condition for zero return curre nt, namely ,f m = 0, is per- 
fectly satisfied inside the beam (jrionda et al.ll2000al ). Also, 
if the mean poloidal (axial) magnetic fi eld, -Bpz, is super- 
imposed to arrange the helica l field fe.g. lAsada et ah 2002; 
iGabuzda fc Chernetskiil [2003 ) , as well as the plasma wall 
in the vacuum channel (Section 1) being in close proximity 
to the beam, the m ode can be stabilized l|Uhm fe Lampel 
Il98d : lDavidsor][l990h . For example, for the rectangular pro- 
file case, the sufficient condition for instability can be written 
as 

f-><^. <*> 

where b (> r^) is the radial distance from r = to the 
plasma wall, q = U) c b/(2&p), and cu c h = eBo/(^bfn c c) is 
the relativistic cyclotron frequency. It turns out, in equa- 
tion (I50|) . that the aforementioned cases of either f m — 0, 
g 2 3> 1 realized for a larger Bo (or a smaller 7b) or r^/b ~ 1 
lead to the robust stabilization. Another possible explana- 
tion is that because the plasma is likely non-uniform in ac- 
tual environments, the mode can be continuously carried 
away from the unstable (to stable) region. However, the 
absolute growing mode tends to be retained at the region 
where the unstable condition is satisfied, thereby the insta- 
bility can evolve (but note again that in the present con- 
text, the growth factor will be favorably small, as reflected 
in equations l47l and l48[) . Note that a similar mode-preference 
is known to appear in laboratory plasmas establishing non- 
uniformity. 

Along with equation (150 ^ , once the aforementioned nec- 
essary steps are taken to stabilize the convective mode, the 
stability properties of the absolute mode also change. Simi- 
lar to the convective mode, for example, for the rectangular 
profile case, the larger values of r^/b and / m = have a 
stabilizing influence on the absolute mode. The superimpo- 
sition of Bo 7^ yields different effects: shifting the range of 
0(= uj — kgVb) corresponding to instability relative to the 
B — case, and expanding the b andwidth (uj 2 h + Cop) 1 ^ 2 
of the instability in the fi-space <|Davidsonlll990h . It ap- 
pears that these effects are not significant enough to disturb 
the reasoning that the convective mode must be somehow 
strongly suppressed or stabilized in the actual system. 

In conclusion, a full derivation has been given of the 
eigenvalue equation for the kink-type perturbations super- 
imposed on the relativistic electron-positron beam with the 
modified Bennett and rectangular profiles. The dispersion 
relations have been examined for the convective and abso- 
lute modes of the resistive hose instability. It has been shown 
explicitly that the pair production effects lower the betatron 
frequency Cop and expand the beam filament radius 7V The 
following key results are found. 

(i) The reduction of (bp leads to a reduction of the growth 
factor of a convective mode scaled as fa ~ Cup/vb- 

(ii) The expansion of prolongs the resistive decay time 
for perturbed current Td (oc r 2 ; equation I33[l. resulting in a 



reduction of the growth factor of the absolute mode scaled 
as fa ~ l/(« b Td). 

For typical parameters, we found Cup 3> l/Va, implying that 
the convective growth is dominant, although its enormously 
rapid evolution, which ought to disrupt jets, appears to be 
incompatible with the observational facts because of the pos- 
sible interpretation given above. 

The Bennett profile effects also suppress the convective 
growth, but the growth distance is still much shorter than 
the observed length of jets (even if the positron effects are 
involved). In addition, the effects remove, for the absolute 
mode, the resonant divergence of the growth factor, and 
the growth distance is found to be comparable to (or 
longer than) the length-scale of jets. To summarize these 
consequences, preferentially, (i) actual evolution of the 
convective mode, in a low frequency and/or high resistivity 
regime, is ruled out and (ii) evolution of the absolute mode 
is allowed, more safely for the Bennett profile case. Closer 
inspections by observations, laboratory experiments and 
kinetic simulations are awaited. 
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APPENDIX A: DERIVATION OF MODIFIED 
EIGENVALUE EQUATION 

In this appendix, we derive an eigenvalue equation that gov- 
erns equations 11311) and 11371). According to the method o f 
characteristics (IKrall fc Trivelpiecel Il986l; iDavidsonl fl990). 
we begin by integrating equation ([26} from z — — oo to 
z' = z. Neglecting the initial perturbation at z' — — oo, we 
obtain 

Sfb(p,r,e,z,T) = 



cpi 



dpi 



Vb 



x (v'r£ + %^SA,[r'(z'),e'(z'),z',T]. (Al) 

Using v' r — Vb (dr' /dz') and v' e — r'vb (d6' /dz'), and equa- 
tions (|29p and (|30|) . equation (|A1|) can be written as 

„, \ e Pz 8f°(p,r) f z ,fdr' 8 d0' 8 

Sfb{p,r) = | / Az — — + — — 

c p± 8p± J \ dz' or' dz' 86' 



xSA z [r'(z')] exp |i [e'(z') -e]-i— (z -z)Y (A2) 

In terms of the generic function of 8J-[r'(z'), 0'(z'), z'], 
we utilize the chain rule for differentiation: dSF/dz' — 
[(dr'/dz') (d/dr') + (dO'/dz') (8/89') + 8/8z'\ ST. Then, 
we have 



Sfb(p,r) 



e Pz 8f h \p,r) 



cpi 



dpi 



a 'I d j_- n 

^ dz' v h 



xSA 



[r'(z')]exp[i[e'(z')-e] -i£ (*'-*)} 



(A3) 



When we impose 8'(z' — z) — 6 and r'(z' — z) — r, and 



recall [5A z (r%, = _ 
gives 



0, the integration of equation 



S.fb(p,r) = 



e Pz <9/b(p,r) 



cpi 



dpi 



5A z (r) + i— 

Vb 



dz' 



x5A z [r'(z')] exp ji [0\z') - 6] - (z - z) | 



(A4) 



Substituting equation (|A4I) into equation (|25[) concomitant 
with the replacement of 8/8t — > 8/8t, we obtain 



1 d r d \ 1 4niuja(r) 
r dr \ dr J r 2 c 2 



5A z (r) 



■■ 47re 2 /3bnbA/p air 



*U w F b (H) 

n 



SA z (r) + —I(n,r, U) 

Vb 



where the orbit integral / is defined as 



(A5) 



I(n,r,U) 



d(j> 
2^ 



d(SA z [r'(z') 



:cxp{i[e'( z ')-e]-i£c} 



(A6) 



Here, £ = z' — z, and <j> is the perpendicular momentum 
phase, conforming to p x + ip y ~ p± (cos (j> + i sin <j>) . 

Concerning the unknown function 5A z [r' (z')] in the in- 
tegrand of equation (|A6|) . it is known that, for |w|rd < 1 
(where Td is the resistive decay time; cf. equations [33] and 
I41[) . SA z (r') — 5A z (r)r' /r provides a reasonably good ap- 
proximation in the beam interio r permeated by the magnetic 
field of Bg oc r l|Davidsonll 19901 ). Physically, such a solution 
reflects a rigid displacement of the magnetic field. Putting 
the approximate expression of 8A z (r') into equation ()A6[) . 
the betatron orbit expression r'(z') exp [i(9'(z')] arises in the 
integrand. As usual, this can be expanded by solving the 
equation of motion of {v h {d 2 /dz' ) + uip[r' (z')]}[x' (z') + 
iy'(z')] = 0. After the manipulation, the orbit can be de- 
scribed as 

r'(z')exp [id'(z'j\ = (1/2) [rexp(ifl) +iu^ (px/jbm e ) 

x exp(i</>)] fexp(— iw^C/«b) ~ exp(iw /3 C/«b)] • (A7) 



As for the integrand of the RHS of equation (|A5[) . we invoke 
the density inversion theorem (cf. Section 3.2), which gives 



™ W F h{ H) = 



dipbjr) 
dr 



1 dn° b (r) 
fib, dr 



(A8) 



Making use of equations (|A7[) and ()A8[) . equation ()A5 



Suppression of hose instability 



(involving equation I A6f) can be integrated, to finally yield 



T he form of equati on (IA9|I is similar to that of equation (45) 



equation (|A9|) we have the corrections related to the pair 
production rate, that is, the factor A/ pa i r and the correction 
via uj% given in equations (|12p or (|22fl . 

This paper has been typeset from a TfjjX/ I^Tj^X file prepared 
by the author. 
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